UNCLASSIFIED 


DEFENSE  DOCUMENTATION  CENTER 

FOR 

SCIENTIFIC  AND  TECHNICAL  INFORMATION 

CAMERON  STATION.  ALEXANDRIA.  VIRGINIA 


UNCLASSIFIED 


NOTICE :  Wien  government  or  other  drawings,  sped- 
fications  or  other  data  are  used  for  any  purpose 
other  than  in  connection  with  a  definitely  related, 
government  procurement  operation,  the  U.  S. 

Government  thereby  Incurs  no  responsibility,  nor  any 
obligation  whatsoever;  and  the  fact  that  the  Govern¬ 
ment  may  have  formulated,  furnished,  or  In  any  way 
supplied  the  said  drawings,  specifications,  or  other 
data  Is  not  to  be  regarded  by  Implication  or  other¬ 
wise  as  In  any  manner  licensing  the  holder  or  any 
other  person  or  corporation,  or  conveying  any  rights 
or  permission  to  manufacture,  use  or  sell  any 
patented  Invention  that  may  In  any  way  be  related 
thereto. 


ig  &93 


D1 -82-0316 


BOE I NG 


SCI  ENTI  FIC 
RESEARCH 
LABOR  AT  OF?  I  ES 


A  Note  on  Regular  Perturbation  Theories 

> 

1 

F.  E.  Bisshopp 


November  1963 


Dl-82-0316 


BOEING  SCIENTIFIC  RESEARCH  LABORATORIES 
FLIGHT  SCIENCES  LABORATORY  REPORT  NO.  84 


A  NOTE  ON  REGULAR  PERTURBATION  THEORIES 


by 

F.  E.  Bisshopp 


November  1963 


A  NOTE  ON  REGULAR  PERTURBATION  THEORIES 


by 

F.  E.  Biashopp 


I.  INTRODUCTION 

There  exist  in  the  literature  several  related  formulations 
of  regular  perturbation  theory*  (c.f.  Morse  A  Feshbach;  Kumar)*  In 
deciding  which  formalism  to  apply  to  a  specific  problem  we  must  take 
into  account  the  radius  of  convergence  of  the  various  methods,  and, 
if  possible,  the  ease  of  application.  In  this  note  we  shall  develop 
a  formalism,  closely  related  to  Fredholm's  resolvent,  which  has  an 
infinite  radius  of  convergence*  The  present  formulation  shares  with 
the  formalism  derived  directly  from  the  Fredholm  resolvent  (Morse  * 
Feshbach)  the  disadvantage  of  being  rather  cumbersome  in  its  applies-* 
tion.  Consequently,  we  shall  go  on  to  develop  an  asymptotic  approxi¬ 
mation  to  the  result*  By  this  means  we  shall  obtain  a  relatively  simple 
procedure  for  evaluating  the  result  of  a  formulation  which  is  unrestricted 
as  regards  radius  of  convergence*  Finally ,  we  shall  discuss  the  close 
relation  between  the  present  theory  and  Feenberg's  (1948)  formulation 
of  regular  perturbation  theory. 
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II.  FORMULATION  OF  THE  CLASS  OF  PROBLEMS 

We  shall  confine  our  attention  to  eigenvalue  problems  which 
may  be  written  in  the  fora 


Vn 


E  t 
n  n 


(1) 


where: 

1.  H(  ia  a  self  adjoint,  linear  operator  which  depends  in 
a  continuous  manner  on  a  complex-valued  parameter  e; 

2*  {E^}  is  a  complex-valued  discrete  spectrum; 

3«  ♦  is  an  element  of  a  Hilbert  space,  Q; 

4-  With  the  possible  exception  of  a  finite  number  of  iso¬ 
lated  singularities  in  the  e-plane,  is  a  bounded 
operator. 

It  should  be  noted  that  the  definition  of  the  space  &  generally 
involves  specification  of  certain  boundary  conditions  on  we  shall 
consider  the  case  where  the  boundary  conditions  may  depend  upon  e,  but 
not  upon  En* 

Now,  the  object  of  a  perturbation  theory  is  to  develop  a 
technique  for  solving  for  *  and  Eq  in  terms  of  known  solutions,  u^  and 
&n,  of  a  related  eigenvalue  problem, 


u  =  &  u 

on  n  n 


(2) 


The  operator  Hq  will  be  assumed  self-adjoint,  and  the  above  problem, 
called  the  unperturbed  problem,  is  to  be  related  to  the  original  in  the 
following  sense: 
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1.  It  must  be  possible  to  define  the  difference! 
Hi  -  Hq  -  Hfi  ,  on  $  ; 


2.  The  space  Q  spanned  by  the  solutions  {u^}  must  be  a 
Hilbert  space  for  which 
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H  ♦  e  $ 
o  o 

HJ  e  $ 
1  o 

E  ♦  e  & 
o 


where  E  is  any  complex  number; 

3«  There  must  be  a  value  of  e  (which  we  may  take  to  be 
e  =  0)  for  which  =  0. 

It  will  be  convenient  to  take  {uq}  to  be  a  normal  set,  i.e* 


(v 


u) 

IB 


am 


(3) 


where  (u,v)  is  a  suitably  defined  inner  product  on  It  should  be 

noted  that  we  have  made  no  specific  assumption  of  independence  of  Hq  on 
e,  though  in  practice  this  is  often  the  case* 

In  virtue  of  the  identity, 

H  =  H  -  (H  -H  )  s  H  -  H,  ,  (4) 

e  o  0  6  oi 


we  may  write 

♦a  '  (Ho  -  En>_1  Vn  •  <5) 

If  we  define  a  projection  operator,  Pn,  according  to  the  property. 


6,  6 

xn  mn 


(6) 
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we  hare 


where 


and 


Note  that  g(e)  merely  fixes  the  normalization  of  t  ;  therefore  it  is 
to  a  certain  extent  a  function  which  may  be  chosen  at  our  convenience* 
It  will  be  convenient  to  take 

g(0)  =  1  (9) 

for  then,  if  the  eigenvalue,  of  the  unperturbed  problem  is  non¬ 

degenerate,  we  have 

K  u  -*  0  and  ♦  -»  u  as  e  -*  0  .  (10) 

n  n  n  n 


III.  BRILLOUIN-WIONER  FORMALISM 

This  perturbation  theory,  perhaps  the  easiest  one  to  derive, 
follows  upon  setting 


g(e)  =  1  and  (1-K)-1=Y(K) 

q  t-L.  n 


*  m«b 

The  eigenvalue  Eq  then  follows  by  computation  of  the  inner  product  of 
Equation  (1)  and  uq;  thus 


E  =  fu  ,  (H  -  H,  )  Y  (K  )m  u  ^ 

n  Vn  0  1  m^b  n  V 


By  definition 


*  («,.  VV.) 


l  n 


where  HXm  5  (V  Hlu«)  * 


where 


H  j  H .  .  ...  H, 


(H  K®  j  =0  for  m  4  0 
°  n/nn 

H  i  H  *  *  • • •  Hj 

vVjiViV 
V“/  (%-En)-**(VE- 

(i)‘  *  i' t -i 

n  X^/n 
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e  =  e  -h 

n  n  nn 


y  /  v  Hni!H  V2  * *  *  *lmn 

vJ  (VE»): "  (VE°) 


The  above  equation  gives  Eq  implicitly;  in  general  the  right-hand  aide 
must  be  truncated  and  Eq  solved  for  numerically* 

The  usefulness  of  this  method  is  often  impaired  by  the  fact 
that  the  domain  of  the  e-plane  within  which  it  converges  may  be  too 
small  to  include  the  problem  under  consideration*  In  virtue  of  the 


identity. 


-'■&  ■ 


.1  -  ^  ■ 


it  follows  that  a  sufficient  condition  for  convergence  is 


*  Max 


((kbu,  Knu>y 


where  the  maximum  is  that  obtained  as  u  varies  over  the  set  {u  3*  Thus, 
if  there  is  a  simply  connected  domain  of  the  e-plane  which  contains  e  =  O 
and  within  which  ||  ||  <1,  the  Brillouin-Wigner  result  will  converge 

within  that  domain  regardless  of  the  details  of  the  dependence  of  Hg 
on  e.  In  practice  it  is  difficult,  if  not  impossible,  to  determine  the 
domain  of  convergence  (even  when  is  uncomplicated)  because  of  the 

appearance  of  EQ(e)  in  Kn* 

It  should  be  noted  that  if  &n  is  a  degenerate  eigenvalue  hav¬ 
ing  associated  solutions  uq,  v^,  ••*,  wq  and  if  the  limit  as  e  -»  0  of 
one  or  more  of  the  quantities 
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(u 


Hnw  ) 

In 


e  -  e 

n  n 


is  nonzero^  then  the  present  formalism  is  inapplicable*  Since  the 
generalization  which  includes  such  cases  is  well  known,  we  shall  not 
treat  such  degenerate  perturbations  in  this  or  any  of  the  following 


sections 


IV.  THE  FREDHOLM  RESOLVENT 


In  this  section  we  shall  discuss  a  formalism  which  is  similar 
to,  but  not  identical  with,  that  based  on  Fredholm's  resolvent  (cf. 

Morse  4  Feshbach,  p.  1018  ft.).  It  will  simplify  matters  considerably 
if  in  what  is  to  follow  we  employ  a  device  by  means  of  which  the  case 
where  depends  upon  e  in  an  arbitrary  manner  is  reduced  to  the  case 

where  H^  is  a  linear  function  of  a  new  parameter,  X.  First,  let  us 
reconsider  what  we  want  out  of  our  perturbation  theory:  the  point  is, 
we  want  a  continuation  of  the  solution  of  Equation  (l)  which  is  given 
by  the  Brillouin-Wigner  theory  near  e  =  0  to  a  region  of  the  e-plane 
which  includes  the  'final'  value,  e^  (say).  In  general  such  a  con¬ 
tinuation  will  depend  upon  the  path  in  the  e-plane  along  which  the 
solution  is  followed;  i.e*,  there  will  in  general  be  branch  points  and 
singularities  of  H  •  However,  once  a  path  0eQ  is  decided  upon  along 
which  is  bounded  and  varies  continuously,  we  may  consider  at  each  e 
lying  on  0tQ  the  related  problem  where 

H  -  X  Hj  .  (18) 

The  process  of  continuation  of  the  solution  along  0eQ  may  now  be  looked 
upon  as  a  continuation  of  the  solution  along  a  path  01  in  the  X-plane 
at  each  value  of  e  on  0eQ.  Since  the  perturbation,  XH^,  varies  linearly 
in  the  X-plane,  the  solution,  is  an  entire  function  of  X,  and  the 

continuation  is  independent  of  the  particular  path  01.  Furthermore, 
if  e  be  taken  sufficiently  close  to  zero,  the  application  of  the 
Brillouin-Wigner  theory  to  the  related  problem  has  a  radius  of  conver¬ 
gence  greater  than  1.  Thus,  we  obtain  a  unique  continuation  along  0cq, 
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satisfying  the  requirement  of  reduction  to  the  Brillouin-Wigner  result 
for  sufficiently  small  e. 

Let  us  now  turn  to  the  continuation  of  solutions  in  the  X- 
plane.  Equation  (7)  becomes 

*„  -  <1-XK„)-1 "»  (19) 

where  we  have  again  set  the  function  g  which  fixes  the  normalization  of 

♦  equal  to  one* 
n 

The  treatment  of  Equation  (19)  given  by  Morse  &  Feshbach 
parallels  the  treatment  of  integral  equations  of  the  Fredholm  type 
(c.f.  Mikhlin,  Section  9)*  The  essential  notion  involved  in  this  treat¬ 
ment  is  that  of  introducing  a  function  x(X)  whose  reciprocal  represents 
the  singularities,  which  according  to  Fredholm's  alternative  appear  at 
values  of  X  for  which  the  homogeneous  equation, 

(1-XK  )«p  =  O  ,  (20) 

m  n  m 

has  a  nontrivial  solution.  The  determination  of  X(X)  is  carried  out 
by  writing  Equation  (19)  as 

*n  =  X^T  [x(X)(1-XKn)“1]  un  ,  (21) 

where  it  i6  required  that  x(X)(l-XK^)  1  shall  be  an  entire  function  of 

X  and  ||  x(X)(l  -  XK  )_1  ||  /  0.  It  follows  that  X(X)  is  an  entire  func- 
n 

tion  of  X  which  vanishes  only  at  the  points  X^  for  which  Equation  (20) 

holds;  the  order  of  the  zero  at  X  is  given  by  the  degree  of  degeneracy 

m 

of  X  • 
m 
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We  shall  not  pursue  this  line  of  development  any  further 
since  detailed  accounts  are  available  in  the  references  cited  in  this 
section*  The  point  we  wish  to  make  is  that  a  somewhat  simpler  for¬ 
malism  results  if  we  observe  that  the  function  X(X)  need  not  be  intro¬ 
duced  in  the  present  problem  since  we  have  at  our  disposal  the  normal¬ 
izing  function  g(e).  We  shall  take  the  normalizing  function  as  g(e,X) 
with  g(e)  =  g(e,l),  and,  suppressing  for  the  moment  the  dependence  of 
g  on  e,  write 

*n  =  g(X)(l  -  XlO”1  un  (22) 

with  the  requirement  that  g( X) (l  -  XKq) shall  be  an  entire  function  of 
X  and  that  ||  g(X)(l  -Xl^)"1  ||  /  0. 

Note  that  the  present  formulation  implies  that  if  ♦  is  to 
be  nonsingular,  then  its  projection  onto  uq  vanishes  at  the  eigenvalues, 
X  .  Thus  it  may  be  seen  that  the  singularities  encountered  in  previous 
treatments  follow  from  the  unrealistic  assumption  that  +n( X)  is  never 
orthogonal  to  +^(0). 

The  function  g(X)  is  not  by  any  means  uniquely  determined  by 
Equation  (22)  and  the  associated  restriction;  a  given  solution  may  be 
multiplied  by  any  nonvanishing,  entire  function  of  X*  The  formal  ex¬ 
pression  for  4n  of  course  depends  explicitly  on  the  choice  of  g(X); 
different  functions,  g(X),  lead  to  different  weighting  of  the  various 
orders  in  the  ordering  of  terms  in  the  solution  according  to  powers  of 
X*  Consequently,  the  number  of  regular  perturbation  theories  we  may 
derive  in  this  manner  is  limitless* 


11 


Now  suppose  we  specify  the  normalization  Integral, 

<tQ,  *n>  =  f(X)  (23) 

According  to  the  properties  of  g(X),  f(X)  may  be  taken  as  any  non¬ 
vanishing,  entire  function.  The  most  natural  choice  for  f(X)  is  of 
course  f(X)  =  1;  however,  we  shall  not  adopt  any  specific  choice  of 
normalization.  Thus  we  have 


[g(X)]2^(l  -  uq,  (l-XKn)_1un^  =  f(X) 

and 

e.  ■  (*»• 

(d-MiJ-1  ...  „J 

(<l-’u£.>‘1  v  0 


(24) 


(25) 


It  should  be  noted  that  when  E^  is  expressed  as  in  Equation  (25),  the 
quantities  f(X)  and  g(X)  do  not  appear,  i.e.,  all  of  the  different  per¬ 
turbation  theories  derivable  by  choice  of  g(X)  (or  f(X))  give  precisely 
the  same  expression  for  E^.  Now  the  above  equation  gives  Eq  as  an  entire 
function  of  X;  therefore  if  we  consider  sufficiently  small  values  of  X 
and  use 

<30 

(1-XK  r1  =  Y  (XK  )m  for  |X|  <  X  (26) 

n  4-  n  o 

m=0 


The  result,  if  expressed  as  a  power  series  in  X,  holds  for  all  values 
of  X.  Thus  we  obtain 
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-  £'^B/£'*fflim  . 


m=0 


m=0 


(27) 


where 


m 


-  S'  <v  K  H„K'  »„> 

1=0 

m-1 

-  S'  (un*  K+iH1Km-1-t  un) 
1=0 
m 

=  y '  (u  ,  k+V~*  u  ) 

k-i  n’  n 

i=0 


(28) 


To  obtain  the  power  series  for  E  ,  we  now  treat  (2xml  )  1  as  the 

n  in 

m 

generating  function  of  a  series  G^.  From  Equation  (13)  it  follows 

that 


( K+m)  =  (K®)  .  0 

nn  nn 


for  m  /  0 


(29) 


and  thus  that 


1=1 

o 


h  "  ° 


(30) 


The  recursion  relation  for  the  quantities,  G^,  is  now  found  to  be 
m-1 

Gjj  =  S  IkGm-k-l  for  m  >  2 

k=l 


(31) 


with  G  =  1,  G^  *  0.  The  power  series. 


e  s  y  \m  F 

n  m 

fflsO 
m 


y  J,  G  , 

t—i  k  m-k 

k=0 


(32) 


with 
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is  then  valid  for  all  values  of  X. 

Equation  (32)  evaluated  at  X  =  1  now  gives  En(«)  implicitly 
at  any  point  in  the  complex  e-plane  which  may  be  connected  to  the 
origin  by  a  curve  along  which  is  a  bounded  operator.  The  drawback 
involved  in  the  use  of  this  formula  is  that  at  X  =  1  it  is  in  general 
necessary  to  retain  a  large  number  of  terms  in  order  to  obtain  a 
reasonable  approximation  to  Eq.  In  the  next  section  we  shall  discuss 
an  asymptotic  approximation  to  the  characteristic  equation  ^Equation  (32)y. 
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V.  ASYMPTOTIC  APPROXIMATION  OF  THE  CHARACTERISTIC  EQUATION 

Let  us  reconsider  Equations  (27)  and  (28)*  What  we  shall 
do  is  to  find  asymptotic  approximations  of  both  the  numerator  and  the 
denominator  of  the  characteristic  equation. 


The  reason  for  keeping  both  numerator  and  denominator  of  the  above 
equation  is  that  when  the  numerator  formally  diverges  the  denominator 
diverges  in  the  same  manner*  As  we  shall  see,  it  is  possible  to  write 
the  asymptotic  approximations  in  such  a  manner  that  the  divergences 
cancel  and  a  well  defined  characteristic  equation  results. 

Now  from  Equations  (15)  and  (28)  we  have: 


for  mil,  where  the  vectors  u  and  u  and  the  matrix  L  are  defined  as 


-  15  - 


n*. 


u 


AXn 


Vv\ 


'/Jt* 


*,  *’  *  n 


(35) 


(in  the  above  definitions  the  branch  of  \^-Eq  is  to  be  assigned  in 
any  consistent  manner,  e*g> 


0  s  ang 


(Vv\)<  = 


(36) 


As  will  be  seen  presently  the  results  are  independent  of  Equation  (36)< 
Also, 


m+1 


£Yk+*H  k”+1-A  -  EYK+iH.Km-^ 

fal  \  °  / nn  A=0\  1  / 

H  ,  ...  H,  ft.  t,  \ 

(  n-t,  l  n  I  l-  *  \ 

. VV 


H  •  •  •  •  Hf  n 


U  .  •  •  •  U|  u 

(ab  nX_  i 

2  )  t - r - 7— 2 — c 

1  \li  n)  U  V 


E  I  , 
n  m+1 


/  +Tm-1  \ 
[u  L  ul 


(37) 


for  a  ^  1# 

The  use  of  Equations  (35)  and  (37)  together  with 


i  =  1  .  In  =  0  ,  j  =  e  f  J-  .  -H 


o  n 


nn 


(58) 


16  - 


gives  the  characteristic  equation  the  form 


(E  -ft  +  XH  )  +  X2  E  Xa(u+L®u) 
n  n  nn _ m=0 _ 

1-X2  Exm  E  aV,  %  'y(u+Lmu) 
m=0  jVn  a(W 


=  0 


(39) 


(It  is  interesting  to  note  that  setting  the  numerator  of  Equation  (39) 
equal  to  zero  gives  exactly  the  result  which  would  be  obtained  by  the 
Brillouin-Wigner  theory. ) 

Now  Equation  (39).  as  it  stands,  still  represents  an  order¬ 
ing  of  terms  in  the  characteristic  equation  in  powers  of  X.  Since  we 
wish  to  take  X  =  1  this  is  clearly  the  wrong  ordering,  and  consequently 
we  propose  to  order  terms  in  Equation  (39)  according  to  quantities 
whose  smallness  is  independent  of  the  value  of  X.  The  small  quantities 
we  shall  use  are  in  fact  the  elements  of  u+,u  and  L.  In  order  to  make 
the  point  clearer,  let  us  make  the  transformation 

x^it)  1  <1,0) 


where 


X 

o 


Max  „ 
l  ,£•  alV 


(41) 


(One,  but  not  both,  of  l  and  l'  may  be  n. ) 

Equation  (39)  is  formally  unchanged  by  the  above  transformation,  but 
the  matrix  elements  and  H^,  are  now  bounded  above  by  1  and  the 

original  problem  is  recovered  for  X  =  Xq.  Now  for  finite  values  of  Er 
there  will  be  an  integer,  l Q,  for  which 


lv«i  -*»  MV*. 

i  i+1 


when  l.  >  t 
i  n 


(42) 
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Thus  we  may  attempt  to  formulate  an  approximation  based  on  the  smallness 
of  the  matrix  elements  of  L  which  follows  from 


"j IV 


for  X*  i  <  >  tn 


(43) 


In  what  follows  we  shall  speak  only  of  the  case  where  n  =  0; 
the  generalization  to  cases  where  n  /  0  involves  essentially  nothing 

new* 


Then 


Now,  suppose  for  the  moment  that 

|WEo|  ‘  IVv^ol 


(v*.  ♦  ™oo>  -  °M] 

1  J0  (x%>m}  * 0 


=  O 


(44) 


(45) 


provides  an  asymptotic  approximation  (valid  in  the  limit  of  large  separ¬ 
ation  of  unperturbed  eigenvalues)  to  the  characteristic  equation-  In 
order  to  evaluate  the  result  of  the  asymptotic  approximation  at  \  =  Xq, 
we  find  the  analytic  continuation  of  both  the  numerator  and  the  denom¬ 
inator-  In  the  case  at  hand  this  is  a  particularly  simple  task  since 
we  need  only  sum  the  series, 


Z  (*%)“ 

m=0 


7^7 


(46) 
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Equation  (45)  now  becomes 

=  0  (47) 

From  the  above  equation  (47)  it  is  clear  that  the  pole  at  X^  such  that 
Eq  =  &o  -  is  not  a  singularity  of  the  characteristic  equation  even 

though  X^  is  (within  our  approximation)  a  point  at  which  the  Brillouin- 
Wigner  theory  diverges.  If  we  now  multiply  through  by  one  troublesome 
factor  and  forget  about  the  denominator  we  obtain  the  simpler  character1 
istic  equation, 

(Eo  -  eo  *  X  Boo)  (Eo  -  ei  *  «ll)  -  *  “('T^)]  <43> 

(Note  that  within  the  accuracy  implied  by  the  above  equation  there  is 
no  value  X^  for  which  the  Brillouin-Wigner  formula  diverges  unless 
^ol^lO  =  75110  apparent  anomaly  is  peculiar  to  the  lowest  approxima¬ 

tion  of  this  type.) 

Equation  (48)  represents  only  the  lowest  order  of  approxima¬ 
tion  in  this  formulation.  If  it  should  turn  out  to  be  necessary  to 
evaluate  the  characteristic  equation  at  a  value  of  X  for  which  Equation 
(44)  is  violated,  we  would  naturally  wish  to  retain  more  terms  in  our 
asymptotic  approximation.  It  turns  out  to  be  possible  to  provide  the 
analytic  continuations  necessary  for  the  evaluation  of  the  asymptotic 
approximations  of  arbitrarily  high  order.  Consider,  for  example,  the 
case  where  terms  up  to  and  including  0  f( cn-Eq )”1"j  are  retained;  then 


20  - 


21 


-  22  - 


It  can  be  shown  that  an  equivalent  way  of  writing  the  1*  approxima¬ 
tion  is, 

„  3°h\°  A-  a0i1\i2%°  \ 

• '  tf. i  “7"“  *  i^i  . 

1  Jl  1  2  Jl  3lJ2  I 

",  a-  H  •"  V 

+  ***  .A,  “HO  Tn5  \  (6] 

3'ml  ’r1*  -%-h  / 


V**jp 


.  ^  ^pjp+i  ^p+i^p 

> '  -3»r- 

P+1  J1"*JP+1 


subject  to  the  restrictions  jm  ^  j^  when  m  >  1.  The  above  form  is 
exactly  the  result  obtained  if  Feenberg’s  approximation  is  truncated 

by  deleting  every  term  which  contains  an  index  greater  than  N;  in  the 
limit  where  N  -*  •  the  two  results  are  the  same. 

Equations  (55)  through  (61)  have  been  derived  by  simply 
neglecting  the  denominator  of  Equation  (52).  Consequently,  they  are 
subject  to  the  appearance  of  apparent  divergence  whenever  AN  -•  0.  How¬ 
ever,  as  discussed  previously,  multiplication  of  the  foregoing  equations 
by  Ajj  to  give  the  result, 

(u^MjjUjj)  =  0  when  Ajj  =  0  ,  (62) 

is  a  justifiable  procedure.  Thus  we  have  a  means  of  justifying  the 
use  of  Feenberg’s  formula  beyond  the  points  where  it  diverges  (where 
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A  «  0,  c*f*  Morse  4  Feshbach  on  this  point)  and,  if  desired,  a  pro- 
n 

cedure  which  may  be  used  in  the  neighborhood  of  such  points  as  well* 

The  convergence  of  our  asymptotic  approximations  is  of  course 
not  insured  in  general*  A  discussion  of  sufficient  conditions  for 
convergence  would  involve  specification  of  the  behavior  of  the  unper¬ 
turbed  eigenvalues,  the  degeneracies,  wq,  and  the  matrix  elements, 

H  ,  as  m  and  n  become  large*  Such  a  discussion  is  beyond  the  scope 

All  “ 


of  this  note 
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VI.  AN  EXAMPLE 

As  an  example  of  the  manner  in  which  the  formulae  of  the  pre¬ 
ceding  section  may  be  applied,  let  us  consider  the  Mathieu  equation, 

2 

M  -  ecos2x  t  ■  Et  ,  HO)  =  ♦( 2tt)  (63) 

dx 


We  take  as  the  unperturbed  problem, 


.2 

— T  =  6u  ,  u(0)  =  u(2n)  , 

dx 


(64) 


with  even  solutions, 


/2rf 

_1_ 

✓w 


cos  nx 


(65) 


and  add  solutions, 


/W 


sin  nrtx 


C  =  -  n 

n 


(66) 


2 

Since  (=  cos  x)  is  symmetric  about  x  =  tt  the  odd  eigenfunctions  do 
not  affect  the  perturbation  of  the  even  ones  and  vice  versa.  Let  us 
fonfine  our  attention  to  even  eigenfunctions;  then  we  need  consider 
only  a  nondegenerate  problem  with 
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From 

the 

formulae 

of  the 

last 

section 

and 

the  ’zeroth’ 

appro  ximat ion 

<Mo  * 

0) 

we  obtain  the  following  results: 

e 

-  E(0> 

-  E(l) 

- 

B  (exact) 

0 

o 

o 

o 

0 

0 

0 

0 

0 

4 

2 

1.5505 

1.54487 

L. 54486 

36 

18 

~  5-8 

~  6 

100 

50 

~  10.8 

~  10 

The  interesting  point  is  that  the  2nd  approximation,  which  in  this  case 
gives  the  relatively  simple  cubic  characteristic  equation. 
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gives  a  result  which  at  e  =  4  is  two  orders  of  magnitude  more  accurate 
than  the  corresponding  quartic  approximation  derived  from  the  Fredholm 
resolvent  (cf.  Morse  4  Feshbach).  In  addition,  at  e  =  100,  a  value 
20  times  greater  than  the  radius  of  convergence  of  the  Brillouin-Wigner 
formulation  of  this  problem,  the  approximate  result  is  still  within  20$ 
of  the  exact*  In  fact,  as  has  been  noted  by  Morse  4  Feshbach,  Feenberg's 
formula  gives  an  exact  continued  fraction  representation  of  the  eigen¬ 
value  in  this  case. 


It  might  be  argued  that  we  have  not  given  the  method  a  fair 
test  in  taking  a  problem  for  which  it  converges.  This,  however,  would 
miss  the  point  —  the  usefulness  of  this  method  lies  in  the  ease  of 


application  on  the  one  hand  and  the  estimate  of  the  error 


("('».! -0)’ 


on  the  other.  Needless  to  say,  the  error  estimates  can  be  refined 
considerably  in  any  specific  problem- 

Finally  it  should  perhaps  be  noted  that  the  computation  of 
our  second  approximation  requires  the  evaluation  of  at  most  nine  dis¬ 
tinct  matrix  elements.  This  is  to  be  contrasted  with  the  evaluation 
(in  principle)  of  infinite  series  required  in  the  more  cumbersome  theory 


based  on  Fredholm's  resolvent 
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